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Abstract 

In the context of the recent interest in solvable models of scattering mediated by 
non-Hermitian Hamiltonians (cf. H. F. Jones, Phys. Rev. D 76, 125003 (2007)) 
we show that the well known variability of the ad hoc choice of the metric which 
defines the physical Hilbert space of states can help us to clarify several apparent 
paradoxes. We argue that with a suitable 9 a fully plausible physical picture of 
the scattering can be recovered. Quantitatively, our new recipe is illustrated on an 
exactly solvable toy model. 
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1 Introduction 



Whenever one considers the one-dimensional differential Schrodinger equation 

ip{x) = E ip^x) , X G (— cxD, cxd) (1) 

in the scattering regime, i.e., with the boundary conditions describing the transmitted 
and reflected waves, 

f e'^"" + C 6-'''% X < -1, 
^(x)= . ' ' (2) 

[ D e'^"" , X > 1 

one usually assumes that the flow of probability is conserved, ICp+lDp = 1. 
Recently, Jones pointed out that several serious conceptual difficulties can arise 
when certain current tacit assumptions (demanding, typically, the reality of the 
potential V{x)) are, tentatively, weakened. He also presented several persuasive 
arguments why one should try to weaken these assumptions. 

Among the latter arguments, the most persuasive support of theoretical as well 
as conceptual innovations would certainly lie in the undeniable success of the re- 
cent, phenomenologically motivated transition to certain manifestly non-Hermitian 
Hamiltonians generating real and observable spectra of bound states [2J. The idea 
(well known to mathematicians [3j) found, recently, several interesting applications 
in nuclear physics [1] and in field theory [5]. 

Of course, there exists an obvious difference between the bound-state problem 
(for which wave functions ip{x) are localized) and the scattering scenario (where all 
the waves remain non-negligible all along the whole real axis). Jones [1] even came 
to a rather sceptical conclusion that the preservation of a sensible probabilistic in- 
terpretation of a generic non-Hermitian model of scattering may be quite costly and 
difficult even when the tentative introduction of a suitable non-Hermiticity in the 
Hamiltonian itself remains restricted to a very small domain of x. Similar observa- 
tions have been also made in a few older scattering models where the violation of 
the rule |Cp + = 1 has been explained and interpreted as a phenomenologically 
acceptable manifestation of the presence of some "hidden" degrees of freedom in the 
model [6]. 
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dx' 



+ V{x) 



Being not satisfied by these "effective" theories the author of ref. [T] formulated a 
much more ambitious project where the physical picture of scattering would parallel 
the above-mentioned "fundamental" theory of bound states based on non-Hermitian 
Hamiltonians [H [3] . Unfortunately, the quantitative results of ref. [Ij were not too 
encouraging (cf. also their recent completion [7]). In essence, the sensible proba- 
bilistic interpretation of the models under consideration seemed to require that the 
above-mentioned standard boundary conditions should be modified to read 

^(x)= : (3) 

[ D e*''^' + D' e"*''^ , x > 1 . 

Unfortunately, this formula contains a strongly counterintuitive "backwards-running" 
component proportional io D' ^ Q in the scattered solution. 

In what follows we intend to weaken the resulting skepticism. We shall start from 
the same basic theoretical premise and postulate that the effective theories of refs. [6] 
are in fact not too interesting since they just mimic the presence of certain unknown 
dynamical mechanisms admitting, i.a., an annihilation and/or creation during the 
scattering. In this sense, we intend to search now for a new quantitative support 
for the possible feasibility and consistency of the alternative fundamental approach 
where one tries to re-establish the conservation of the probability. 

In section [2] we shall commence our considerations by a brief and compact review 
of the most relevant results of ref. [1] . We summarize there the overall philosophy of 
the fundamental theory where the input Hamiltonian H [say, of eq. ([T])] is interpreted 
as a mere auxiliary operator. One assumes that with this auxiliary operator the 
formal calculations become exceptionally simple. At the same time, the "correct" 
physics is assumed to be defined, via an invertible map f2, in terms of a certain 
"true" physical Hamiltonian 

h = nHn-^ (4) 

which is expected to be prohibitively complicated. For illustration one can recollect 
the "true nuclear physics" of ref. |1| where h was a full-fledged Hamiltonian of an 
atomic nucleus while H has been constructed as its much more easily tractable map. 
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In all the similar scenarios, the similarity mapping Q must be assumed non- 
unitary - otherwise, one would return to the mere traditional, Hermitian class of 
Hamiltonians. In order to suppress or circumvent the related, mostly purely technical 
obstacles we decided to employ our recent bound-state experience [5] and to restrict 
our attention to a fairly restricted class of the solvable illustrative dynamical models. 
They are introduced and described in section [3l 

The manifestly non-perturbative character of our present class of models will 
enable us to make use of the flexibility of the mappings Q and to draw several con- 
sequences from the exact solvability of our models. In section H] we shall show how 
some of the apparently unavoidable paradoxes of ref. [I] can find their explana- 
tion and resolution. In particular, for our manifestly non-Hermitian set of specific 
Hamiltonian-representing operators H ^ W we shall demonstrate that they can be 
assigned a consistent and unitary physical interpretation of the scattering based on 
standard asymptotic boundary conditions ([2]). 

In section |5] a compact summary of our present constructive arguments will be 
complemented by a few optimistic remarks concerning the possible future extension 
of the class of our present models of scattering towards some less schematic non- 
Hermitian Hamiltonians. 

2 Scattering from complex potentials 

The mathematical background of the fundamental models of scattering from localized 
non-Hermitian centers will be illustrated here on a set of solvable models. On this 
level we shall demonstrate that a very natural interpretation of this type of scattering 
is feasible. On an abstract phenomenological level we shall stress that in our present 
update of the extension of the scattering theory of ref. P] a core of observational 
consistency should and can be sought in an at least partial, asymptotic survival of 
the observability of the coordinates. 
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2.1 Jones' solvable example 

For the majority of the real and smooth one-dimensional short-range potentials 
the description of the scattering is routine. One solves the ordinary differential 
Schrodinger equation under the standard asymptotic boundary conditions. For the 
complex V{x)s in ([1]), fundamental theory makes the scattering unitary via an ap- 
propriate adaptation of the metric G in the Hilbert space of states. Unfortunately, 
an unexpected and unpleasant consequence has been detected in pLj where the re- 
placement of the standard asymptotic boundary conditions ([2]) by their fairly coun- 
terintuitive "amendment" ([3]) has been found necessary in principle. Fortunately, in 
a concrete illustration using the complex delta function toy potential of ref . [H] , 

yiMostafazadeh) = 2A(1 + Is) 6{x) (5) 

the contribution of the nonvanishing coefficient D' to the flow of probabilities proved 
negligible [1]. Thus, in the leading-order approximation it was possible to return to 
the original boundary conditions ([2]). 

After such an approximate confirmation of the internal consistency of the funda- 
mental theory approach a careful perturbation analysis of scattering by potential ([5]) 
has been performed in [1] leading to the intermediate result 

In the naive effective-theory interpretation of this formula the non-conservation of 
the fiow of probability merely refiects the fact that the manifestly non-Hermitian 
Hamiltonian H is merely an auxiliary operator defined in the "wrong" Hilbert space 
^(unphysicai) _ Yilis to employ the map Vt to move to another, unitarily non- 

equivalent "correct" Hilbert space 7-^(p^J/s«caO ^y^q^q the "true" representant (jlj) of the 
Hamiltonian remains safely Hermitian [1]. 
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2.2 Long-range non-localities induced by the short-range 
potential: a paradox 

One of the most unpleasant formal features of physical metric 6 = Q{H) in the 
Hilbert space of states T-^Cp'iy^^caZ) i\ig^t it is usually strongly non-local even if the 
original potential is local, V = V{x) (cf. also [8J). This implies that it is hardly 
feasible to perform any computations in "^iphysicai) ^ ^j^g purely auxiliary Hilbert 
space 

y^iunphysicai) ^j^g computations are assumed much easier. All the ket vectors lib) 
which lie in the latter space lack, unfortunately, any direct physical interpretation. 
Even the standard requirement of the observability of the Hamiltonian degenerates, 
in this auxiliary space, to the identity [1] 

H^ = QHQ-K (7) 

This relation represents the Hermiticity of the Hamiltonian h in '}{ipf^y^^<^°-'') ^ 

h = QHQ-^ = h^ = (Q-^y WQK 

We may deduce that one has to put = Q'^Q [11]. In the language of mathematics 
we must guarantee that the physical metric 6 is compatible with relations ([7j). Vice 
versa, any operator a = cr^ representing an observable in "^ivhy steal) j^^^g have its 
appropriate quasi- Hermit ian partner S in 'j-^i'^^physicai) -^j-^idi obeys an analogue of 
eq. ([7]) using the same metric [12] . 

The shift e in eq. has been assumed small in ref . [I] . This opened the possibility 
of an explicit use of Mostafazadeh's metric [9] available in perturbation-series form 

Q(Mostafazadeh) ^ ^ _ J ^ ^ ^(1) ^ Q(^^2^ _ ^g) 

In coordinate representation the unit operator / becomes represented by the delta- 
function kernel 6{x — y) but a manifest and large nonlocality emerges already in the 
first perturbation order. 



ry«(x,?/) = ^sgn(?/2-a;2) [^(xy) e^^l^"^! + ^(-xy) e^'l^+^'IJ . (9) 

The emergence of such a nonlocality leads to serious problems because "the physical 
picture of the scattering is completely changed" since, on the positive half- axis the 
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physical wave function "no longer represents a pure outgoing wave . . . but . . . contains 
an 0{6) component of an incoming wave as well" [1]. 

We shall be able to show that while the deformations caused by the use of the 
locally non-Hermitian interaction remain long-ranged in their character, they need 
not necessarily lead to the emergence of spurious components in the outcoming wave. 
Such an observation is not in contradiction with the fact that "one should change 
the Hilbert space by adopting the appropriate metric [which] must differ from the 
standard one not only in the vicinity of the non-Hermitian potentials, but also at 
distances remote from it" [Ij. Nevertheless, we shall argue that at least some of 
the spuriosities emerge only due to an inappropriate choice of a specific metric, the 
definition of which is known to contain infinitely many free parameters [ll|TTl|T3]. In 
this sense we shall make use of a simpler model and recommend here the construction 
and use of another, quasi-local metric operator O = O^'^-^) ^ rj. 

3 Discrete Schrodinger equations 

Some of the standard scattering-theory considerations can be simplified when one 
replaces the ordinary differential equation ([1]) by the difference equation 

ji^ hV{xk)^p{xk) = Eip{xk) . (10) 

For example, in some pragmatic numerical calculations one chooses a sufficiently 
small step-size h > and introduces discrete coordinates, 

Xk = kh, k = 0,±l,.... (11) 

This makes the usual real line replaced or approximated by an infinitely long discrete 
lattice. The most elementary application of such a discretization occurs when one 
wants to construct bound states. For certain real as well as complex potentials a 
sample of the construction may be found in our papers [Tl] . Some of them illustrate 
also the fundamental-theory approach to the non-Hermitian quantum bound states 
where the Hamiltonian H is treated as quasi- Hermitian, i.e., Hermitian only in the 
Hamiltonian-dependent Hilbert space 7^(p'^s/««ca0 
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3.1 Discrete in and out free waves 



Let us assume that the potential in eq. ffTOj) vanishes beyond certain not too large 
distance from the origin, V{x±j) = 0, j = M, M + 1, . . .. In the free-motion domain, 
we abbreviate ipj = ^iid 2 cosip = 2 — h'^E and replace eq. (fTOj) by recurrences 



ifjf}, + 2 cosypz^f - = 



or by the matrix equations Hoip^'^^ = h'^k'^^^^^ or Mq(lp)iP^^^ = 0, viz., 
/ -. •. •. : \ / : \ 



(12) 



'•. 2cosv9 —1 
— 1 2cos(/9 —1 
—1 2cosv9 



V 



do) 



(13) 



which may be assigned a doublet of independent solutions. 



= const ■ g^, Q± = exp(±i (p) . 

Precisely in the spirit of ref. ^ we can speak about a PT— symmetric free Hamilto- 
nian. 



2 -1 
-1 2 -1 
-1 2 



V 



It is easy to verify that its spectrum is real, provided only that the new energy 
variable ip is real [8]. This imposes an inessential constraint —2 < 2 — h'^E < 2 
upon the energy range, i.e., we must have E G (0,4//i^). At any finite choice of 
the lattice step h > this is reminiscent of the similar feature of the spectra in 
relativistic quantum systems. This connection has been given a more quantitative 
interpretation in ref. [T5] . 
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Let us finally add an interaction with nonzero elements forming merely a finite- 
dimensional submatrix in the Hamiltonian. The scattering of an incoming wave may 
be then characterized, say, by the boundary conditions 



e'""^ + Re- 



imip 



m < -M , 
m > M - 1 



i.e., by eq. ([2]) in its discrete version. 



(14) 



3.2 Discrete short-range model of scattering 

In a way inspired by our recent studies of finite-dimensional non-Hermitian Hamil- 
tonian matrices H with real spectra [HI [IS] let us pick up one of these models and 
contemplate its infinite-dimensional generalization which would admit scattering so- 
lutions. Its explicit matrix representation will be tridiagonal, one-parametric and 
doubly infinite. 



2 -1 
-1 2 
-1 



Hi 



(15) 



-1 + a 2 -1 

-1 2 -1 
-1 2 



Inside the interval of a G (—1, 1) all the 2 /T— dimensional truncations of Hi have 
the 2K— plets of eigenvalues which are all real and lie inside the closed interval 
(1,3). In the non-Hermitian regime with the growing \a\ these energies pairwise 
degenerate at the "exceptional points" a = ±1 and get complex at \a\ > 1. At K = 2 
and K = 3 the smoothness of the a— dependence of these truncation-dependent 
standing wave energies has been illustrated in [8]. At K = 1 their explicit form 



reads h'^k^ = 3 ± — a'^ and offers a schematic guidance and a nice quantitative 
illustration of what happens in general. 



In a search for the transmission and reflection amphtudes our infinite-dimensional 
matrix problem for scattering Hi ip^^^ = h'^k'^ tp^^^ can be split in its two free-motion 
parts ffT^ for the respective "in" and "out" solutions (fUl) valid up to M = 1. They 
have to be matched near the origin, 



-1 e'^ + e-'^ -1-a 
-1 + a e'^ + e-'^ -1 



T 



1 




1 — a 



-Ivp 



These two linear equations for R and T can be simplified, 



1 




-(1 + a) 






R + 




T = 


1 - a)e'^ 









-1 



It is easy to write down their explicit solution. 



;i-a)(l-e2i^) 



A 



A = l-(l-a2)e' 



This gives an exact analogue 

l-a[l + f/(a,(^)]-i 



(16) 



l + a[l + f/(a,(/?)]-i ' " ^""-^ 2(l-a)(l-cos2(/?) 

of equation Nr. (11) of ref. pj. In both cases, the sum of probabilities is greater 
than 1 or less than 1 depending on the sign of the deviation of the coupling constant 
form its Hermitian zero limit. The same conclusion can be read in ref. [1] so that 
in the weak-coupling regime our present difference-operator parameter a plays the 
same dynamical role as its differential-operator predecessor e in eq. ([5]). Moreover, 
due to the nonperturbative character of our result one can rewrite eq. f|T6|) in the 
equivalent form [17] 

(17) 



|^2M |y|2_ + 4(1 - a)2 sinV 



a4 + 4(1 - a2) sin^ ^ 

which is more compact and clarifies the nature of the singularity reached in the limit 
a ^ 0. 

We shall also see below (cf. section that after the necessary adaptation of the 
Hilbert space of states and after the ad hoc modification of the inner product the net 
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result of the changes will be the elementary modification of the coefficient of sin^ ({> 
in the numerator to 4(1 — a^), thus restoring the usual and physically consistent 
unitarity of the scattering. 

3.3 More-parametric models 

With another Hamiltonian 



2 -1 






-1 2 


-1-6 




-1 + 6 


2 -1 - a 






-1 +a 2 


-1-6 




-1 + 6 


2 -1 






-1 2 ■-. 



one can expect that many results obtained previously for its 6 = special case H\ 
can find a natural generalization. The same expectations concern also the next-step 
candidate with three free parameters, 





-1 2 -1-c 
-1 + c 2 


-1-6 






-1+6 


2 

-1 + a 


— \ — a 
2 


-1-6 






-1 + 6 


2 -1-c 
-1 + c 2 -1 



etc. In all of them the asymptotic "in" and "out" solutions of eq. f|T^ remain 
uninfiuenced by the interaction. Equally well, the matching of these in and out 
solutions remains feasible at any number k of parameters. 
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0, 



In our first nontrivial scattering model H2 ip^'^'^ = h'^k'^ ip^'^'' the matching may be 
mediated by the choice of M = 2. This means that the following four matching 
conditions must be considered, 

e-^'^ + Re'^'^ 

e-'^ + Re'^ -x-i 
T + xo 

-1 Te'^ 

The first line defines the correction X-i and the last line defines the correction xo, 

{l + b)x-i = b (e-"^ + R e'n , {l-b)xo = bT. 

This reduces the number of our equations to two again, 

e-^''^ + Re^'^ 

-1+6 2cosy9 — 1 — a 
— 1 + a 2cos(p —1 — 6 



-1 2cos(/? -1-6 

-1 + 6 2cosv9 -1-a 

-1 + a 2cosv3 -1-6 

-1 +6 2cos(/? 



r'^ + Re'^j /{1 + b) 
T/(l-6) 



0. 



After their simplification we may easily eliminate 

(1 + 6)T _ 1 + i?e^'^ 
(1 - a) (1 - 6) ~ TTlP~e^ ' 

We end up quickly with the explicit definition of R for our second model if 2? 

+ 2 62 cos 2 + 6^ 



R 



A 



A = 1 - (1 - a^) e^'^ + 262 ^2i^ ^ ^4:^ _ 



We observe a close parallelism with the preceding model. From the easy first-order 
estimates 

(1 + 6)T 



a) (1 - 6) 



1 + 0{a') + 0{b' 



R = 0{a^) + C(62) , 

we may immediately deduce that 

|i?|2 + |r|2 = l-2a-Ab + 0{a^) + 0{b^). (18) 

This two-parametric dependence parallels closely the one-parametric prediction of- 
fered by eq. f|T6l) . 
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4 Unitarily non-equivalent Hilbert spaces 



In a climax of our paper we shall make use of the fact that our models are exactly 
solvable, at least in terms of the methods based on the computer-assisted symbolic 
manipulations and extrapolations. This will enable us to construct many eligible 
candidates O for the metric in the physical Hilbert space. In contrast, even the 
construction of their subclass denoted by the symbol p and possessing an explicit 
perturbation form 

g = I + ^er]^'^ + 0{e^). (19) 

was an achievement for differential operators in refs. [1] and [9]. In such a context, 
a core of our present message is that due to the simplified, difference-operator rep- 
resentation of observables we shall be able to select a better metric Q^'-^^^ ^ rj. With 
its use, some of the most counterintuitive manifestations of the nonlocality paradox 
will simply disappear. 

4.1 Ambiguity problem 

In the effective interaction scenario, formulae ([6]) or (fT6l) and (ITSl) would certainly 
indicate the presence of an absorption and/or creation at e 7^ or a 7^ and b ^ 0. 
In the fundamental theory one assumes a change of the Hilbert space such that the 
original (i.e., standard) definition of the inner product 

= [ ^*{x)^'{x)dx (20) 
Jm 

is replaced by its more general weighted version in the new space, 

(M')q= [ ip*(x)e(x,x')'^'(x')dxdx' . (21) 

The purpose of such a change is in making the Hamiltonian selfadjoint. 

It is well known that the choice of the inner product (^T^ is ambiguous [4J. One 
of the standard constructive solutions of the ambiguity problem giving = 1] has 
been proposed by Mostafazadeh [9]. In the mathematically most easily tractable 
dynamical regime of a very small deviation \6\ <^ 1 from Hermiticity this author 
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arrived at the explicit perturbation approximation ([8]) + ([9]) where the maxima of 
function ri{x, y) lie on the two perpendicular lines defined by the trivial equations 
X ± y = in the x — y plane. Subsequently, the latter recipe has been used in 
ref. PP where operator = g was defined as a self-adjoint square root flTI?]) 

of metric t] emphasizing that in terms of physics, "the relevant wave function is not 
^(x) = {x\4)), but ^(a;) = (a;|*) = {x\g\^y. 

In this context, the mathematical essence of our present amendment of scattering 
theory lies precisely in an innovation of the choice of Q and O since among all the 
available mappings Q the selected g remains also very strongly non-local, indeed. 

4.2 The existence of diagonal matrices G = G^*^^^ 

In the technically most complicated part of our present considerations we decided to 
choose a Hamiltonian and to treat its quasi-Hermiticity condition (J?!) as a linear set 
of equations for all the matrix elements of the metric. 

In the first attempt we choose Hi and verified that there exists the infinite- 
dimensional matrix solution Gi of eq. ([7]) which is diagonal, i.e., in our present 
terminology, quasi-local. 



1 — a 




1 - a 




1 — a 






1 + a 




1 + a 




1 + a 



This result was obtained via tedious symbolic manipulations on the computer. Its 
simplicity is both very surprising and very useful because one of the integrations in 
the related inner product fl?I]) drops out. Moreover, its diagonal kernel can trivially 
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be factorized into the product of two diagonal operators p = v9, i.e., 



SQL) 
Pi 



vT+a 



VI + a 



They remain self-adjoint and positive definite at all the not too large real as. 

The diagonality of the latter matrix enables us to insert it in eq. Nr. (17) of 
ref. [1] and to deduce that the explicit formula for the "correct" operator X of the 
observable coordinate coincides with its standard diagonal-matrix form with elements 
given by eq. ( ITTi) above. In the same manner one can also recall eq. (jll) and introduce 
the operator h^^^^ = p^*^^^ Hi (^p^^^^^ which represents the isospectral Hermitian 
Hamiltonian of our system and which replaces eq. (fTSll by the real and symmetric 
tridiagonal matrix 



2 -1 
-1 2 -1 
-1 2 




-VI -a2 


2 -1 






-1 2 -1 






-1 2 ■ 











f22) 



We see that this operator differs from the purely kinetic Hamiltonian just in a small 
vicinity of the scattering center. In spite of such a strict locality of the interaction, 
the metric itself remains deformed far away from the scattering center. 

The fact that the manifest non-Hermiticity of our toy model Hi did not involve 
the mixing of incoming and outgoing waves that occurred in the model of ref. [1] 
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encouraged us to proceed towards the more complicated models using the same 
brute-force method. After the choice of the next, two-parametric non-Hermitian 
Hamiltonian H2 the calculations still remained sufficiently easy for us to deduce 
and verify the existence of the following two-parametric quasi-local solution Q'^^^ of 
eq. ([7]) represented by the diagonal matrix 



(1 - a) (1 - bf 




(1- a) (1-6)2 




(1 - a) (1 - b') 






(1 + a) (1-62) 




(1 + a) (1 + 6)2 




(1 + a) (1 + 6)2 



Similarly, we took A; = 3 in the next continuation of the series of solutions 0^ 
pertaining to Hj^. It is easy to verify that the three-parametric quasi-local solution 
e^^^^ of eq. dZD is still obtainable diagonal matrix with the same elements 
(1 — a) (1 — 6)2 (1 — c) 2 in all the upper left corner and with the similar array of 
the same elements (1 + a) (1 + 6)2 (1 + c)2 in its lower right corner. The remaining 
"central" quadruplet of the "anomalous" diagonal elements is formed by the following 
four-dimensional diagonal central sub matrix of our doubly infinite matrix 03*^^'*, 

" (l-a)(l-6)2(l-c2) 

(l-a)(l-62)(l-c2) 

(l + a)(l-62)(l-c2) 

(l + a)(l + 6)2(l-c2) _ 

The general pattern of extrapolation is now obvious. It would be easy to write down 
and, via eq. ([7]), verify an immediate extrapolation of the k = 1, k = 2 and k = 3 
matrices G^'^^'* to the higher subscripts k whenever necessary. 

There are all reasons why, in the context of scattering, the specific diagonal 
metrics B^,*^^^ should be preferred in comparison with all their non-diagonal and, 
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hence, more nonlocal alternatives. With such a new postulate we may now return 
to paper [1] once more. First of all our results reconfirm the high plausibility of 
the hypothesis that even the violation of the Hermiticity which is strictly localized 
in space should be expected to influence, manifestly, even the asymptotics of the 
wave functions. The explicit analysis of our schematic models indicates that even 
our minimally nonlocal metric operators remain, strictly speaking, different from the 
most common Dirac's delta-function metric Qi^^^"-'^) (^x , y) = 5{x — y) at all distances. 

This being said, we found it quite fortunate that at the sufficiently large distances, 
i.e., for |x| ^ 1 and/or \y\ ^ 1, the difference between 6(^*™^)(a;, and Q^^^\x,y) 
degenerated, in all of our models, to the mere introduction of a nontrivial multipli- 
cation factor, 

Q^Q^\x,y) = const{signx) ■ e^^"''"'\x,y) , |x| > 1 , > 1 . (23) 

In another formulation, our explicit constructions very strongly support the affirma- 
tive answer to the question of the existence of a "spatially localized non-Hermiticity" . 
A formal key to such an answer is that in a schematic model we constructed certain 
new and very specific, "quasi-local" metrics B*-*^^^ with the property fl23|) . 

The strict validity of this proportionality rule at almost all the coordinates x 
and y may be admitted to be an artifact resulting from our specific tridiagonal- 
matrix choice of our "toy" Hamiltonians iJfc. Still, the validity of such a rule at all 
the sufficiently large coordinates may be expected to survive transition to a larger 
family of models and, perhaps, also to some slightly less friendly generalized quasi- 
linear forms of B*^'^'^'', say, with a band-matrix structure. This would still allow us to 
conjecture that with the metrics B = B*-'^^-* the internal consistency of the models of 
scattering (and, in particular, of their asymptotic boundary conditions) would not 
be violated after the extension of the present theory towards many less schematic 
and reasonably non-Hermitian models of dynamics. 

One of the instrumental versions of our conjecture will have the form of the 
requirement = in boundary conditions so that eq. ([3]) = eq. ([2]). Then, the 
choice of a unique metric B*^'^'^-' characterized by its minimalized nonlocality should 
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be perceived as strongly recommended in the conceptually consistent fundamental 
scattering theory using non-Hermitian Hamiltonians. 

Marginally, the latter requirement can be supported also by the remark that 
"it has been known for some time that ... for the potential ix^ and the infinite 
PT— symmetric square well ... the particle is confined ... so that the range of the 
non-locality is limited. Scattering potentials highlight this feature [of nonlocality] to 
its full extent because the wave functions ... do not have compact support." [Ij. This 
means that the "traditional" choices of 7^ QiQ^) can still offer a fully consistent 
model of the physical reality for bound states. After all, we already noticed that 
many models with O 7^ 0(Q^) found applications in nuclear physics [4j and in field 
theory [5J. Other constructions of G 7^ QiQ^) appeared also in the coupled-channel 
problems [12] or in the Klein- Gordon- type models [18] etc. 

5 Summary 

In the differential-equation model of scattering studied in ref. [T] the behavior of 
the physical "in" and "out" states was strongly non-local so that, for example, the 
outgoing waves contained a non-negligible "incoming" component. In this context 
we showed here that such a paradox is not inevitable for non-Hermitian systems 
with real spectra. A set of counter-examples has been described here in which the 
local non-Hermiticites carried by the Hamiltonian implied just the necessity of the 
replacement of the usual scalar product (l20l) by the local, re-scaling change of the 
measure, 



This enabled us to address several conceptual difficulties as encountered in ref. [T] 
where the description of scattering caused by several short-range non-Hermitian sam- 
ple potentials V{x) has been presented. In this context we discovered and described a 
family of non-Hermitian short-range Hamiltonians Hi, H2, ... for which the descrip- 
tion of the scattering looks almost as easy and natural as in the standard Hermitian 
regime. 




(24) 
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Our selection of short-range interaction models proved technically much simpler 
than expected. We revealed several amazingly close parallels with their continuous 
delta-function analogues. We were able to bring new arguments, first of all, thanks 
to certain "unreasonable efficiency" of our non-perturbative method. In this frame- 
work, our main mathematical result is that the hermitizing metrics = 0^ which 
we attached to Hi, H2 and H3 and, by an easy extrapolation, to any Hk are all 
represented by the (infinite-dimensional) diagonal matrices. We believe that this is 
not just a friendly feature of our specific models but rather a generic property of the 
metrics since one has a lot of freedom of their modification in general. 

In the latter spirit, our present main recommendation is that for all the realistic 
non-Hermitian models of scattering one should still try to insist on the requirement 
that the physical metric B is not too non-local. In our present text we succeeded in 
supporting the latter recommendation by a series of the explicit illustrations of its 
feasibility. One of reasons was that we choose the discretization of the real line of 
coordinates as our principal methodical tool. 

The first hints offering a background for such a decision were already found and 
formulated in [S] . The present results can briefly be characterized as a successful tran- 
sition from the bound-state models (or, formally, from the flnite A^— point lattices 
of ref. [Ml) to the scattering scenario (or, formally, to the limit N 00), comple- 
mented by the replacement of the simplest possible one-parametric model of ref. [8] 
by the whole set of dynamically nontrivial localized Hamiltonians Hj^ containing, in 
principle, an arbitrary flnite number k of coupling constants. 

A formal beneflt of our choice of the models appeared to lie in their two-faced solv- 
ability. Its flrst face was rather technical and concerned an easiness of construction of 
the reflection and transition coefficients. Certain massive cancellations in the linear 
algebraic matching conditions made the final formulae unbelievably compact. The 
second friendly face of the solvability emerged during our systematic construction 
of the metrics O. An easiness of the guesswork encountered during extrapolations 
k —>■ k + 1 is worth mentioning since it proved helpful and saved computer time. 

A priori we couldn't have hoped in the amazing diagonality of our solutions 
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of eq. ([7]) or in their asymptotically constant form or in a "user-friendliness" of 
the transition from the trivial model Hk with = 1 to virtually all of its > 1 
descendants. We firmly believe that at least some of these properties will also be 
encountered in some other, similar but less schematic models of the dynamics. 

Needless to add that many emerging questions remain open. Some of them (like, 
typically, the numerical efficiency of the discretizations and an analysis of the practi- 
cal rate of their convergence) have been skipped intentionally. The omission of some 
other points was only made with regret, mainly because of their lack of any imme- 
diate relevance for physics. For example, a marginal but interesting benefit of the 
discretization with a fixed gap h > could have been seen in the emergence of paral- 
lelism between the discrete-lattice formulae and their continuous-limit counterparts. 
Besides such a direct possible correspondence between HkS and point interactions, 
another correspondence (viz., to the truncated, finite lattices) has also been omitted 
as too mathematical, in spite of its potential relevance for the verification of the 
reality of the spectra. 

We are sure that even within the domain of physics we did not list all the open 
questions. Pars pro toto, let us sample, in the conclusion, the possible relevance of 
the present models with the extremely simple metrics in the context of the path- 
integral formulation of Quantum Theory where an extremely interesting discussion 
just appears in print [19], concerning the questions of the role of the explicit form 
of the metric G in the partition functions Z[J] and in certain related formulae in 
thermodynamics and/or quantum field theory. 
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